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1. Introduction 

Understanding the quantum chromodynamics (QCD) phase diagram is one of the most inter- 
esting problems in quark and hadron physics. Since the lattice QCD Monte-Carlo (MC) simulation 
has a sign problem at finite chemical potential }x [|]], it is necessary to invoke some approximations 
in QCD or to apply effective models in order to elucidate the whole phase boundary, including the 
position of the critical point The strong-coupling lattice QCD (SC-LQCD) is one of the most 
instructive approximations especially to investigate the phase structure at finite temperature T and 
chemical potential ji. Based on the successes of SC-LQCD in pure Yang-Mills theory [^], chiral 
transition at finite T and ji has been investigated at strong coupling [Q, ^, ^, For example, the 
phase diagram structure has been predicted in the strong coupling limit (SCL) [^, ^, and it is re- 
cently confirmed qualitatively in MC simulations ^ based on the monomer-dimer-polymer (MDP) 
formalism [^]. In order to make a step forward towards the true phase diagram, it is necessary to 
develop the formalism to include the plaquette effects both in MC simulations and SC-LQCD. 

In this proceedings, we evaluate the effective potential in SC-LQCD with one species of un- 
rooted staggered fermion at finite T and ji. including the next-to-leading (NLO, ^(1/g^)) [[TOl, [TT[] 
and the next-to-next-to-leading (NNLO, i^(l /g^)) [ [l^ ] effects, and investigate the coupling depen- 
dence of the phase diagram and the critical point. 



2. Effective Potential in NNLO SC-LQCD 



The effective potential ^gff is obtained in two steps in a finite T treatment of SC-LQCD. 
We first integrate out the spatial links Uj and obtain the effective action ^eff. In the next step we 
integrate out the quark field x and the temporal hnks Uo, and obtain the effective potential ^gff- 
We define the effective action and potential {S^s and ^eff) on a lattice with spatial (temporal) size 
L{Nx) at chemical potential ji as. 



exp[-5eff(;t,^,t/o)] = j &Uj e-^LQCD = J ^y. ^-s^ 
^[X,X,Uo] exp[-5eff] , 



exp 



2N, 



(2.1) 
(2.2) 
(2.3) 

(2.4) 



where d = 3, mo and rjj^x = (— 1)^""' ' are the spatial dimension, the bare quark mass and the 
staggered phase factor, respectively, 5lqcd = Sf + Sg is the lattice QCD action, and Up denotes 
the trace of a plaquette P. Mesonic composites are defined as = XxXx> = Xx^^^q^Xxj^^ ^"^^ 



yx=x 



x+O" 



0,x 



Xx- 



In order to perform the strong coupling expansion at finite T systematically, the cumulant (or 
coupled cluster) expansion is indispensable. By using the cumulant expansion, the expectation 
value of exp[— ^g] is found to be, 

1 



-SsCL 



mj e'^" e-^° = exp 



1 "1" 



,-SscL 



,-Sf 



(2.5) 
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Figure 1: Diagrams contributing to the SCL, NLO and NNLO effective action terms. Quarks (anti-quarks) 
are shown in open (filled) circles, spatial (temporal) link variables are represented by solid (dotted) arrows, 
and plaquettes are represented by squares. Note that for NLO and NNLO diagrams, we also consider the 
hermite conjugate of these diagrams. 



The bracket (• • • )c is called a cumulant and shows the connected diagram contributions, e.g. (5^)^ = 
{S^^ — (So) ■ We find that the sum in the exponent corresponds to the strong coupling expansion 
of the effective action, 

Ses =SscL - £ (5^),, = 5scL + A5nlo + AWo + ^(1//, l/Vd) . (2.6) 

The «-th term in the sum is proportional to and we can identify n = I and n = 2 terms as 

NLO and NNLO corrections (A^nlo and A^nnlo), and Sscl shows the SCL effective action. 

We consider zero, one, and two-connected plaquette configurations, and integrate out Uj after 
putting fermionic elements (xUx and Z^^X)- the leading order of the l/d expansion we 
find that the diagrams shown in Fig. [I] contribute to the effective action, 

5eff = - V7) - ^ I [MM]j, (2.7) 

.V ^" xJX) 

\k\^j,\i\^j,\i\m 

^'^ x,j>0,\k\^j 



+ ^^^^ , I ^ ^ [V+V- + ([MM] [MM] . ^^j^sj , (2.9) 



We have introduced a short-hand notation, [AB\j^x =AxB^^j, [ABCD]jk,x = AxB^_^_jC^^j^iD^^i. Ad- 
ditional mesonic composites appears as W+ = e^f" XxUo.^Uq x+oXx+26 ^.i = ^^^^Xx+io^lx+o^L^''- 
The coefficients are defined as, = d/2Nc , = d{\ + \/2g^)/N^g^ , jS, = d{d - i)(l + 
l/2g2)/8A^4^2 ^ ^ d/2N^g' , j8,, = d{d-\){d- 2)/l6Njg^ , /3„ = d{d - l)/2N^g^. 

We can convert the above effective action in a spatially local and bi-linear form of quark fields 
by using the extended Hubbard-Stratonovich (EHS) transformation [ |l0| , [TTI ]. The fermion part of 
the transformed effective action is found to be in the form of, 

= ^E(CV+-CV-)+£mM,-j8,,£(vA„W,+ + v/,,W,-) , (2.10) 
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where i/Att and va^t are the auxiUary fields and C, C, m contain auxiliary field contributions. We 
further reduce this effective action by introducing a gluonic dressed fermion, 

X!, = x,-2przVrTe''Uo,,X,+d/C , X'. = Xx + 2^,xWra,+o^'^uljC . (2.11) 

With x' and x' , we can absorb the next-to-nearest neighbor (NNN) interaction terms via up to 
/g^) in the modification of mass and the coefficients of V^, m, C and C. 
We now find that the NLO and NNLO corrections lead to the coefficient modification of y+, 
V~ and M in the fermionic effective action. These modifications are interpreted as the modification 
of the wave function renormalization factor Z^^, quark mass niq and chemical potential p.. The 
effective potential ^eft is obtained in a similar way to that in SCL, and is found to be, 

^eff =^iff^ + r,(m,;/i, r) -NAog ^/Z^ , (2.12) 

111 1 

-^eS =2^'o + 2^^' ^^^-^ + ^^''^^'^ P-crV-crV-cr + Pss%,W'< + l^PrsVtsVrs , (2-13) 
mq=^==, m' = b'^a + mo-prt{¥tz + ¥rt) , fl = n - log y^Z^JzI , (2.14) 

Eq{mq) = arcsinh {niq) , (2.15) 

jSr' =i8r + j8„V^r. , K = Av + ^^ssWss , b'^ = ba + 2A>. + ^l^ssYss + 2j8^,VA„ , (2.16) 

Z+ =1 + j8/v/t + 4j8„m'v/TT , Z_ = 1 + jS/i/A^ + 4j3„mV^T . (2.17) 

The auxiliary fields introduced during the bosonization procedure have to satisfy the stationary 
condition, d^^s/d^ = 0, where represents one of the auxiliary fields. The stationary condition 
of ^eff tells us that auxiliary fields other than a and ft)^ are expressed as functions of (cr, (0^). Thus 
the effective potential is a function ofT , }JL, and two order parameters; the chiral condensate a and 
G)t. We can regard (O^ is a vector potential field for quarks; the chemical potential shift is mainly 
determined by (O^, and co^ contributes repulsively to the effective potential in equilibrium. This 
two order parameter feature may be a natural consequence from the potential term from quarks, 
'fq{mq;jl,T). There are two independent derivatives, dYq/dmq and dYq/dfi, which appear in the 
equilibrium condition, then we have two degrees of freedom. 



riog 



sinh[(jV, + l)£^/r] 
sinhteo/r] 



+ 2cosh(A^cAi/r) 



3. Chiral Phase Transition in NNLO SC-LQCD 

The effective potential determines the vacuum and the phase structure of QCD matter. The 
vacuum is determined by solving the stationary condition of =^eff> which is equivalent to searching 



for the saddle point of ^eft in the (a, Wt) plane In Fig. g, we show the the critical temperature 
at /I = {Tc^p,=q) and the critical chemical potential at T = (/ic,r=o)- In both NLO and NNLO, the 
phase transition on the T-axis is found to be the second order, and Tc^^=q decreases as j3 increases. 
NLO and NNLO give almost the same values of Tc^^=o. We also show the MC results on Tc in 
SCL and the critical coupling (j8c) for given values of A^t = 1/7" at = (filled triangles) pj| ]. 
The suppression of T^-^^=q in NLO and NNLO is not enough to explain the MC results. The phase 
transition on /x-axis is numerically found to be the first-order in NLO in the coupling range studied 
here (j8 < 6), and for j3 < 5.5 in NNLO. The first order phase transition is determined by the 
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Figure 2: The j3 dependences of a critical temperature at /i = and a critical chemical potential at T = 0. 
We also show the MC results of Tc in the strong coupling limit and the critical coupling j3e for a given A^t- 
MC results are taken from |l3]l. 



magnitude relation between Eq and p.. As P increases, both of Eq and p. are suppressed, and the 
suppression effects on p. is slightly larger than that in Eq. As a result, IJ^^^jIq is a slightly increasing 
function of j3 . 

In Fig. ^, we show the j8 dependence of the phase diagram in NLO and NNLO. While Tc^^=o 
and /if,r=o are very similar in NLO and NNLO, the shape of the phase boundary and the critical 
point are sensitive to the NNLO corrections. In NLO, the critical point (CP) temperature {Tqp) 
decreases and the CP chemical potential (jUcp) stays almost the same as j8 increases. The tri- 
critical point at strong coupling starts to deviate from the second order phase boundary at j3 ~ 4.5, 
and it becomes the critical end point (CEP) even in the chiral limit. In NNLO, both Tqp and Hcp 
decreases as j8 increases, and the CP stays on the second order phase transition boundary. 

The CP behavior in NNLO is considered to be an improvement over NLO. Since one species 
of staggered fermion corresponds to Nf = 4, the phase transition at /i = is expected to be the first 
order in the continuum region due to anomaly contributions [14]. In the strong couping region, 
the chiral SU/(4) symmetry is not manifest with staggered fermions. The lattice QCD action 
with staggered fermion is invariant under the chiral transformation, Xx — ^ ^w{i^^x)Xx^ where £^- = 
(^_iyo+xi^ is a 75 related phase [Q]. This transformation corresponds to the 75 (g) 75 rotation in 
the spinor-flavor space in the continuum region, which is only a small subgroup of SUa(4). As a 
result, effective number of flavors is small at strong coupling, and the second order phase transition 
at /X = emerges. If we can take account of the finite coupling effects correctly, the order of the 
phase transition at /i = should change from the second to the first order at a certain coupling 
strength. The behavior of CP in NNLO suggests that the first order transition boundary comes 
closer to the T-axis, and agrees with the above expectation while the shift is not enough. 

Another interesting point is the existence of the region where ij.c^^^\t) < /xi^"'^''(r) is satis- 
fied. Between these critical chemical potentials, /ic'**''(r) < n < /Xc^"'^'(r), chiral condensate is 
suppressed compared to its vacuum value but the chiral symmetry is not fully restored. This par- 
tially chiral restored (PCR) matter suggested in NLO SC-LQCD [10] persists to exist also with 
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Figure 3: Phase diagram evolution in NLO and NNLO SC-LQCD. Phase boundary are shown by thin Hne 
surface, and thick surfaces show the first order phase transition boundary. 

NNLO corrections, and would correspond to the quarkyonic matter suggested at large Nc |[r^]. 
4. Summary 

In this proceedings, we have derived an analytic expression of the effective potential in the 
strong coupling lattice QCD (SC-LQCD) with one species of unrooted staggered fermion for color 
SU(A'c) including the next-to-leading (NLO) and the next-to-next-to-leading (NNLO) effects at fi- 
nite temperature and density, and investigated finite coupling effects on the chiral phase transition 
and the phase diagram. Effective action terms have been systematically evaluated based on the 
strong coupling expansion in the leading order of the \/d expansion. NLO and NNLO effective 
action terms have been found to come from one plaquette and two connected plaquette configu- 
rations. We have applied the extended Hubbard-Stratonovich transformation to bosonize fermion 
interaction terms. The next-to-nearest neighbor interaction appears in NNLO corrections, and is 
evaluated by introducing a gluonic dressed fermion. We have obtained the effective potential as a 
function of temperature {T), chemical potential {pi), and the two order parameters: the chiral con- 
densate ((7) and the vector potential (cOt)- The vacuum is determined from the stationary condition 
of the effective potential with respect to each auxiliary field. NLO and NNLO effects result in 
modification of the wave function renormalization factor, quark mass and chemical potential. 

The critical temperature at /i = {Tc,j^=o) and chemical potential at T = (/ic,r=o) are found to 
be similar in NLO and NNLO. The critical temperature at = {T^ .^=o) is a decreasing function of 
j8 in both NLO and NNLO, while the suppression is not enough to explain the Monte-Carlo results 
at N-c = l/T = 2,4 and 8. The critical chemical potential at T = ilJ.cj=o) is found to stay almost 
constant in both NLO and NNLO, due to the cancellation of the effects from mass and effective 
chemical potential reduction. The shape of the phase boundary and the position of the critical point 
(CP) are found to be sensitive to the NNLO effects. As j3 increases, while CP moves towards the 
;U-axis in NLO, CP moves towards the T-axis in NNLO. 

Inclusion of the Polyakov loop and higher order terms in the l/d expansion would be interest- 
ing directions of study. 
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